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In a recent Comment |jj on our earlier preprint g] 
proposing a theoretical explanation for the observed 
strong temperature dependence of the low temperature 
resistivity of low density "metallic" Si inversion layers ||] 
Kravchenko et al. have raised a number of issues which 
require careful consideration. In this Reply to their Com- 
ment we discuss the issues raised in ref. Jl[| with respect to 
our work []2| and respond to the specific questions about 
our model |^] brought up by Kravchenko et al. 

It was conceded in ref. jjj that our theoretical results 
H "are qualitatively similar to those observed experi- 
mentally" but questions were raised about our choice of 
parameters for the charged impurity density (Ni) and the 
free carrier density (n s — n c ) with the implication that 
our choice of Ni and n s — n c may be inconsistent with 
the experimentally used samples in ref. [[| . We first point 
out that neither quantity (Ni or n s — n c ) has been ex- 
perimentally measured, and therefore the authors of ref. 
[Q are expressing their opinions based purely on theoret- 
ical speculations rather than making statements of facts 
m r cf. Q. (This distinction was not clearly made in ref. 
jij.) Their remark [fij on our choice of the charged impu- 
rity density Ni is based on a naive misunderstanding of 
our theory whereas our choice of the free car- 

rier density (in fact, this is the definition of the effective 
free carrier density in our model |Q) is based on more 
subtle arguments than indicated by the discussion of ref. 

®- 

As emphasized in our original paper [Ej the parame- 
ter Ni only sets the overall resistivity scale in our the- 
ory in the sense that the resistivity p oc Ni — Ni 
does not affect the calculated temperature dependence 
at all. The actual value of Ni is set in our work by 
demanding overall quantitative agreement between the- 
ory and experiment, and the very low value of iVj (~ 



0.3 x 10 



10, 



reflects the anomalously high peak mo- 



bility (~ 71,000cm 2 /Fs at lOOmK) of the Si-15 sample, 
which has a peak mobility roughly a factor of five to 
seven higher than that in other typical good quality Si 
MOSFETs (~ 10,000 - 20, 000cm 2 /Vs). An indepen- 
dent confirmation of the low value of Ni in the Si-15 
sample comes from the empirical formula connecting Ni 
and the peak mobility (/x m ) derived by Cham and Stern 
@ : Mm = 1250(Ayi0 12 )-°- 79 where \i m is expressed in 




cm /Vs and iVj in units of 10 



12, 



This empirical 



formula leads to an iVj 
of two of the value N ? 



0.5 x 10 10 cm 2 within a factor 
0.3 x 10 10 cm -2 used in our work 



H. The actual discrepancy is even less because the em- 
pirical formula Q| refers to the T = mobility, which for 
Si-15 should be considerably higher than 71, 000c?7i 2 /Vs 
because of the strong temperature dependence of the mo- 
bility observed in ref. . Thus the low value of Ni used 
in our work is necessitated by the anomalously high 
(almost an order of magnitude higher than usual) mo- 
bility of the Si-15 sample. The fact that a low value of 
Ni is needed to obtain quantitative numerical agreement 
between our theory and the experimental Si-15 data Q| 
at high densities (n s 3> n c ), where the standard Drudc- 
Boltzmann transport theory should be eminently valid, 
shows rather decisively the correctness of our choice of Ni 
We emphasize that the low value of Ni used 
2| simply reflects the anomalously high quality of 
the Si-15 sample of ref. ||, and nothing else. 

Having established that our choice of Ni in ref. M is 
quite reasonable we should now mention two aspects of 
our model H which seem to have been overlooked in the 
Comment [If. First, characterizing the charged impurity 
scattering by a single two-dimensional (2D) charge den- 
sity of N with all the impurity centers located randomly 
in a plane placed precisely at the Si-Si02 interface, as 
we do in ref. 0], is surely a highly simplified zeroth or- 
der model for a complicated situation where the charged 
impurity centers will be distributed over some distance 
inside the oxide layer. We use the simple model to keep 
the number of free parameters to a minimum (just one, 
Ni). If we modify the model slightly, for example, dis- 
place the random charged impurity plane some finite dis- 
tance inside the oxide (or consider a three-dimensional 
impurity distribution, as is likely), we could considerably 
increase 2Vj, making it sound "more reasonable". We be- 
lieve that such fine tuning is unwarranted in a zeroth 
order model and should be left for future improvements 
of the theory. Second, our estimated Ni is surely a theo- 
retical lower bound on the possible value of the charged 
impurity density because our theory j^] clearly overesti- 
mates the impurity scattering strength as we neglect the 
modification in scattering due to electron binding to the 
impurities, and assume that the charged impurities scat- 
ter via the screened Coulomb interaction without taking 
into account the binding effect discussed qualitatively in 
ref. ||. Thus, the parameter Ni in our model M should 
be taken as an effective (single) parameter which charac- 
terizes the overall impurity scattering strength in the sys- 
tem, which should not necessarily be precisely the same 
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as the number density of bare Coulomb scatterers in the 
sample. 

Now , we discuss the second point raised in ref. 0] re- 
garding our choice of (n s — n c ) = n e as the "effective" 
free carrier density participating in the "metallic" Drudc- 
Boltzmann transport for n s > n c . This is a rather subtle 
issue because in our model Q n e is, by definition, the 
free carrier density entering the conductivity a — n e e/i 
where /i is the carrier mobility — Drude-Boltzmann the- 
ory allows for an intuitive separation of the conductiv- 
ity into an effective carrier density (n e ) and an effec- 
tive carrier mobility (/x) . Our choice of n e = n s — n c as 
the effective free carrier density leads to excellent qual- 
itative (actually, semi-quantitative) agreement between 
our calculated theoretical results and all the published 
experimental results in Si MOSFETs and p-type GaAs 
heterostructures as far as the temperature dependent re- 
sistivity on the "metallic" side (n s > n c ) is concerned. 
Our theory ^] also produces the experimentally observed 
non-monotonic temperature dependence in the resistiv- 
ity arising from the quantum-classical crossover behavior 
which, we believe, to be playing an important role. We 
note that our choice of the effective free 

carrier density is not only consistent with the impurity 
binding/freeze-out scenario 0||, but also with the per- 
colation model of ref. || . 

The important open issue is, of course, a direct ex- 
perimental measurement of n e , the effective free carrier 
density near threshold (i.e., n s > n c ) at zero magnetic 
field. It has been known for long time @ in MOSFET 
physics that a direct determination of the free carrier 
density near threshold at low temperatures is almost im- 
possibly difficult (particularly if n s < 10 n cm -2 as it is 
in the cases of our interest), requiring a number of simul- 
taneous complex measurements which, to the best of our 
knowledge, have never been attempted in the samples 
showing the so-called 2D metal-insulator-transition (M- 
I-T) phenomena. Since there seems to be confusion or 
misunderstanding on this point jlj we take the liberty of 
quoting from the authoritative review article Q by Ando, 
Fowler, and Stern: "If trapping, band tailing, or deple- 
tion charge are important (as they are near threshold) 
no single experiment can unambiguously give the mobil- 
ity and carrier concentrations The interpretation 

of the measurements is increasingly unreliable at carrier 
concentrations below 10 12 cm~ 2 ." (p. 490 in ref. Q.) 
We believe that a direct measurement of the zero-field, 
low-tempcrature, low-density, near-threshold free carrier 
density through careful capacitance studies (which are 
extremely problematic at low densities and low temper- 
atures) is the only way to decisively settle the question 
of the effective free carrier density participating in the 
"metallic" transport in the Si samples of ref. ||. 

Following the procedure used in ref. || for GaAs sys- 
tems, we can however determine the operational value of 
n e participating in the Si-15 sample. We show in Fig. 1 
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FIG. 1. The calculated zero-temperature conductivity 
o(T = 0) for the Si-15 sample of ref. [3] as a function of 
density n s using the theory of ref. [2]. The dots are the ex- 
trapolated T — experimental conductivities at respective 
n 3 values taken from ref. [3]. The dashed line is the best fit 
to the experimental results whereas the thick (thin) lines are 
the theoretical results using a Dingle temperature To = 0.0K 
(0.5K) in the quantum screening (see ref. [2] for details). 

the extrapolated zero temperature conductivity a(T = 0) 
of the Si-15 sample || plotted as a function of the carrier 
density n s , both for experimental data || and our theo- 
retical results [||. It is clear that er(T = 0) cx (n s — n*) 
in Fig. 1 where n* w n c . We believe that this specific 
functional behavior of the T — conductivity, vanishing 
approximately linearly in (n s — n c ) — n e near thresh- 
old, provides compellingly strong phenomcnological sup- 
port for our model of using n e — n s — n c as the effective 
free carrier density in the Drude-Boltzmann theory. Note 
that in our model a = n e efi with n e = n s — n Cl and the 
calculated /i has very weak n s dependence near threshold 
(see Fig. 1) — this is strong support for our assumption 
that n e — n s — n c provides a measure of the effective 
free carrier density in the theory. We have found [[l0| the 
same behavior in all the published Si data. 

Third, we address the issue of the unpublished finite 
magnetic field (e.g. Hall effect) measurement of the free 
carrier density alluded to in ref. This connection 

made in ref. is, in fact, misleading since the applica- 
tion of an external magnetic field to a 2D system com- 
pletely changes its physics by converting it to a quan- 
tum Hall system whose localization properties are still 
poorly understood. The free carrier density (n e ) en- 
tering our Drude-Boltzmann theory is, by definition, 
the zero-field free carrier density which is not necessarily 
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the same as that measured in finite field quantum Hall- 
type experiments alluded to in ref . JO . This fact is most 
easily demonstrated by pointing out that the finite field 
measurements of the type mentioned in ref. j^j give the 
free carrier density to be n s even deep inside the local- 
ized regime (n s < n c ) where the free carrier density in 
the Drude-Boltzmann theory is obviously zero, n e = in 
the insulating phase. Thus the finite-field carrier density 
measurements give the free carrier density to be n s ev- 
erywhere (both in the metallic and the insulating phase) , 
and cannot therefore have anything to do with the pa- 
rameter n e entering our Drude-Boltzmann theory. While 
understanding these interesting finite field measurements 
is clearly an important theoretical challenge (well beyond 
the scope of our manifestly zero field theory ^ aimed 
exclusively at understanding the temperature dependent 
resistivity on the "metallic" side), we disagree with the 
suggestion in ref. that an understanding of the zero- 
field transport properties, as in ref. must somehow 
depend crucially on a complete theoretical understand- 
ing of the quantum Hall behavior of these systems. (None 
of the proposed theoretical models for 2D M-I-T can ac- 
count for the finite field quantum Hall behavior.) 

Finally, we mention that characterization of our the- 
ory as a "simple classical model" in the concluding sen- 
tence of ref. El is misleading. While quantum interfer- 
ence corrections are neglected in our theory (we argue 
that quantum interference corrections are overwhelmed 
by the screening effects considered in our theory in the 
T > lOOmK temperature range considered in the exper- 
iments ||), the striking temperature dependence Q of 
our theoretical results arises entirely from the classical- 
quantum crossover phenomenon and the quantum screen- 
ing of charged impurity scattering - a classical theory 
would not have any of the effects obtained in our calcu- 
lation 
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